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Abstract 

Using the Pauli-Fierz model of non-relativistic quantum electrody- 
namics, we calculate the binding energy of an electron in the field of 
a nucleus of charge Z and in presence of the quantized radiation field. 
We consider the case of small coupling constant a, but fixed Za and 
ultraviolet cut-off A. We prove that after renormalizing the mass the 
binding energy has, to leading order in a, a finite limit as A goes to 
infinity; i.e., the cut-off can be removed. The expression for the ground 
state energy shift thus obtained agrees with Bethe's formula for small 
values of Za, but shows a different behavior for bigger values. 
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1 Introduction 

Quantum electrodynamics (QED) has proved one of the most successful 
theories in physics. One of its striking features was the explanation of the 
so-called Lamb-shift, the splitting of the energy levels of the hydrogen atom. 
Most of the results, however, are of perturbative nature, and very little is 
known concerning rigorous results starting from a well-defined (Hamiltonian) 
theory. 

A few years ago, Bach, Frohlich and Sigal initiated a rigorous non-pertur- 
bative study of non-relativistic QED. They made a detailed spectral analysis 
of the Pauli-Fierz Hamiltonian, which is a model appropriate for studying the 
low energy properties of QED, and which is the model we use in this paper. 
Among other results, they proved in [BFS] the existence of a ground state 
for small values of the coupling constant a and the ultraviolet (UV) cut-off 
A, which has to be introduced in the interaction term between the electrons 
and the photon field in order to obtain a well-defined Hamiltonian. Recently, 
Griesemer, Lieb and Loss [GLL] succeeded in removing these restrictions on 
the parameters. We refer to [GLL] for an extensive list of references on the 
subject. 

Our goal in this paper is to calculate the binding energy of one electron 
in the field of a nucleus of charge Z. In contrast to the case of a Schrodinger 
operator without coupling to the quantized photon field, the "bare mass" 
appearing in the Hamiltonian is not equal to the physical mass. To obtain 
the binding energy as a function of the physical mass, which can be measured 
in experiments, and hence to compare the calculated binding energy with 
the measured one, one has to substitute the physical mass for the bare mass. 
To do this, one first has to calculate the physical mass as a function of the 
bare mass, the coupling constant a, and the ultraviolet cut-off A. We do this 
via the dispersion relation, i.e., the energy of a free electron as a function of 
the total momentum (= electron momentum + field momentum). 

There are three parameters in the problem: the coupling constant a, the 
UV cut-off A, and the nuclear charge Z. We consider the case when a is 
small, but Za not necessarily small. I.e., we first fix A and Za and consider 
a as a small parameter. 

Our calculations are to leading order in a. We give rigorous error esti- 
mates, but do not focus on optimizing the error terms. In particular, we do 
not have sufficient control on their behavior for large A to be able to prove 
the finiteness of the result after removing the cut-off beyond leading order 
in a. This remains an interesting open problem. Steps in this direction were 
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taken in [LL1], where the behavior of the self energy for large A was studied. 
Bounds on mass renormalization and binding energies were recently given in 
[LL2] . 

We shall rigorously show that after renormalizing the mass the UV cut- 
off can be removed, at least to leading order in a. More precisely, the binding 
energy, as function of the physical mass and the coupling parameter a, has, 
to leading order in a, a finite limit as A — ► oo. We thus obtain an expression 
for the ground state energy shift due to the presence of the radiation field, 
which holds to leading order in a, and for all values of Za. It agrees with 
Bethe's formula [Be] for small values of Za, but shows a different behavior 
for bigger values. We emphasize again that in our formula for the binding 
energy no UV cut-off appears, in contrast to Bethe's formula, who used the 
dipole approximation for it's derivation and neglects the spin of the electron. 
In Bethe's formula a logarithmic dependence on A appears. Bethe argues 
that a physically reasonable value for A is mc 2 , the energy needed for pair 
production of electrons and positrons to take place. By inserting this value 
for A in his formula he obtains excellent agreement with experiments, at 
least for the case Z = 1, i.e., the hydrogen atom. 

Shortly after Bethe's work, it was noted in [KL] that the Lamb shift 
ought to be finite after removing the cut-off A if one does not use the dipole 
approximation (see also [AF]). By means of perturbation theory the Lamb 
shift was calculated, for both a and Za small, in [Gr]. 

Our proofs below essentially show that second order perturbation theory 
is correct. I.e., to obtain results up to order a only the addition of one 
photon is needed. We use and extend previous ideas in [H], where results on 
the leading order of the self energy and the binding energy were obtained. 
(Concerning the question of enhanced binding, see [HVV].) 

The paper is organized as follows. In the next section, we describe the 
precise setting and introduce some notation that will be used throughout 
the paper. In Section 3 we calculate the dispersion relation and thus obtain 
an expression for the physical mass in terms of the bare mass, i.e., for the 
necessary mass renormalization. The binding energy is calculated in Sec- 
tion 4. In Section 5 it is shown that after renormalization of the mass, the 
UV cut-off can be removed, and a finite expression for the binding energy 
to leading order in a is obtained. See Theorem 4 on page 867 below for 
the precise statement. Finally we comment on the Lamb shift of metastable 
excited states in Section 6. 
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2 Setting and Notation 



We now describe the precise setting and introduce useful notations. The 
Hamiltonian under consideration is the Pauli-Fierz Hamiltonian 

H = —— Up + yJaA(x)) 2 + yfea ■ B(x)] + V(x) + H f . (2.1) 

It acts on H = L 2 (R 3 , d 3 x; C 2 ) <g> J^(L 2 (M 3 , d 3 k; C 2 )), where T h denotes the 
bosonic Fock space. Self-adjointness of H on an appropriate domain has 
recently been shown in [Hi]. 

Units are chosen such that h = c = 1. The electron charge is then given 
by e = y/a, with a the fine structure constant. In nature a ~ 1/137, but we 
consider it here as an arbitrary, small number. The positive parameter mo 
is the bare mass of the electron. 

The • product stands for the usual scalar product in M 3 , and a denotes 
the vector of Pauli matrices, acting on the spin variable of the electron part 
of the wave function. The operator p = —iV x is the electron momentum. 
The vector potential in Coulomb gauge is 



A(x) 



1 

2^ 



E 

A=l,2 



s x {k) a x {k)e ik - x + a* x (k)e- % *- x d 6 k . (2.2) 



—ik-x 



\k\<\ yj\k\ 



Here the integration is restricted to momenta \k\ < A, i.e., we make a sharp 
UV cut-off. We could allow for a more general, spherically symmetric cut-off, 
but we restrict ourselves to this case for simplicity. The vectors £\(k) G M 3 
are any two possible orthonormal polarization vectors perpendicular to k. 
The corresponding magnetic field B = curl A reads 



B{x) 



2vr 



kf\E\{k) a x (k)e tk ' x -a* x (k)e 



A=l,2' 



l<A 



—ik-x 



d 3 k 



(2.3) 



The operators a\(k) and a x (k) are the usual creation and annihilation op- 
erators for a photon of momentum k and polarization A. They satisfy the 
commutation relations 



[a u (k),a* x (q)] = 5(k - q)5 uX 



(2.4) 



and 



K(k), a* x (q)} = , [a„(fc), a x (q)} = . (2.5) 

The field energy is given by Hf = c2r(|A;|), where we denote, in general, the 
second quantized version of a function h(k) by 

dT(h(k)) = V f h(k)a* x (k)a x (k)d 3 k . (2.6) 

Atiy^ 3 
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A general wave function can be written as a sequence 

\I/ = {Y>o,Y>i,V2,---} , (2.7) 

with ipn = tp n (x, k±, . . . , k n ) an n-photon wave function. We suppress the 
dependence on the spin of the electron and the polarization of the photons 
for simplicity of notation. The creation parts of the vector fields A and B 
create a photon with wave function G(k)e~ lk ' x and H(k)e~ tk ' x , respectively. 
Here G(k) = (G 1 (k), G 2 (k)) and H(k) = (H 1 (k) , H 2 (jfe)) are vectors of one- 
photon states, given by 

and 

Their vector components will be denoted by Hi and Gi, respectively (i = 
1, 2, 3), and likewise generally for vectors in M 3 . The function 6 denotes the 
Heaviside step function, taking the ultraviolet cutoff A into account. 



3 Dispersion Relation 

We first consider the case of a free electron, i.e., V = 0. The Hamiltonian 
is then translation invariant, i.e., it commutes with the three components 
of the total momentum P = p + dT(k). It is therefore possible to write the 
Hilbert space and the Hamiltonian as a direct integral 



= / d 3 PH P (3.1) 

JR3 



n 

and 

r® 

H= d 3 PH P , (3.2) 

with Hp acting on Hp. Each Hp is isomorphic to C 2 ®^. The Hamiltonian 
Hp in this representation is given by 



H P = 



2rriQ 



(P - dT(k) + y/aA(0)) 2 + \faa -B(0) + H f . (3.3) 



For any fixed total momentum P, we define £p to be the ground state energy 
of Hp. The existence of a corresponding ground state is of no concern to us. 
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In [C] it was shown, at least for the spinless case, that there does not exist a 
ground state of Hp for P / unless one imposes an infrared cut-off on Hp. 

The function £p is often referred to as the dispersion relation. For small 
P, £p ~ £q + \P\ 2 /2m, where m is, by definition, the physical mass. It is 
strictly bigger than the bare mass ijiq. For given A, mo has to be chosen 
such that m is equal to the mass of the electron, which is determined by 
experiment. The bare mass goes to zero as the UV cut-off A goes to infinity, 
which is usually referred to as renormalization of m . 

We now calculate this mass renormalization to leading order in a. We 
restrict ourselves to the case P < tuq, since we are only interested in the 
behavior of £p for small P. For large P the behavior is different due to 
the different energy-momentum relations of the relativistic photon and the 
non-relativistic electron. 

THEOREM 1 (Dispersion Relation). For \P\ < mo, the energy £p of 
a free particle with total momentum P satisfies 



2rriQ 2ir 2 mo J\ k \ < _^_ \k\ 2 + Ifel — P • k/mo 

' ' — 2mQ 

+\P-k\ 2 ( ', m2 - tIt)) + 0(« 3/2 ) • (3.4) 



|fc| 2 1 1 \(\k\ 2 + \k\) 2 \k\ 



Strictly speaking, we do not claim that the error term is truly 0(a 3//2 ), 
but only that it is bounded by const, a 3 / 2 for small a. The same remark 
applies to Theorems 2 and 3 below. 



Proof. Throughout the proof we will set mo = \ for simplicity. The correct 
dependence on mo is easily obtained by dimensional analysis. We will provide 
appropriate upper and lower bounds on £p. 

We start with the upper bound: We choose as a trial state 



T, 



—\/a 



|fc| 2 + |A:| -2P-k + a 



2P ■ G(k) + a ■ H(k) 



1,0,0,... 



(3.5) 



where G(k) and H(k) are defined in (2.8) and (2.9), and where f is a short- 
hand notation for the vector (1,0) in spin space. The additional factor a 
in the denominator of the one-photon part serves as an infrared cut-off to 
make the norm of ^ finite. 



C. HAINZL, R. SEIRINGER 



853 



We write A = A(0) = D + D*, where D denotes the part of A containing 
only annihilation operators. A straightforward calculation, using that 

1 



A 2 = (D + D*) = -A 2 + 2D* ■ D + D* ■ D* + D ■ D 



7T 



(3.6) 



and k ■ D* |= k ■ G{k) |= 0, yields 
{$\H P \V} = (|P| 2 + ^A 2 ) (*|*) 
- a({2P ■ G + a ■ H) | 



(2P • G + a ■ H) T 



\k\ 2 + \k\ -2P-k + a 

+ a 2 ((^\^>) - l) + const, a 2 . (3.7) 

Using the fact the G is real and H is purely imaginary, and the anti- 
commutation relations for a, namely 



(Jidj + GjOi 



25al, 



ijJ-C 2 > 



(3.8) 



we get 

'(2P-G + a-H) T 



\k\ 2 + \k\ -2P-k + a 

E 1 < H 



A=l 



(2P • G + a ■ H) T J = 
A\P -G x {k)\ 2 + \H x {k)\ 2 



\k\ 2 + \k\ -2P-k + a 



(3.9) 



We now insert the expressions (2.8) and (2.9) for G and H and use that 
Za=i \ P ' £ A(fc)| 2 = \P\ 2 - \P ■ k\ 2 /\k\ 2 by the orthogonality relations of k 
and e\. Moreover, since {^\^) = 1 + 0(aln(l/a)), we arrive at the upper 
bound 



Sp < \P\ 2 + -A 2 

7T 



a 



(2vr) 2 



, 3 4(|P| 2 '^ 



|P-fc| 2 ) +2|fc| z 



|fc|<A 



|fc| 3 (|£f + |fc| -2P-A;) 

+ 0(a 2 ln(l/a)) 



(3.10) 



Lower Bound: We first show that, for suitable constants Q > 0, 

#p> (|P| 2 + |dT(A;)| 2 )(l-Cia) + C , 2ir / -C , 3a . (3.11) 

Together with the upper bound (3.10) this shows that, in an approximate 
ground state * , (^o|^/|^o) < const, a and (* | \dT(k)\ 2 |^ ) < const, a. 

The fact that e\{k) ■ k = implies that dT(k) ■ A = A- dT(k). Using this 
we can estimate by Schwarz' inequality 

± V^(ip\(P - dT(k)) ■ A\il>) = ±2^^{i/j\(P - dT(k)) ■ D\if>) 

< aa{ip\\P - dT(k)\ 2 \ip) + ^\D*D\^) (3.12) 
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for any positive constant a. This leads to the operator inequality 



±^(P - dT{k)) ■ A < aa\P - dT{k)\ 2 + ^-D*D . (3.13) 

In the same manner, we obtain 

±y/aa ■ B < ^E*E + aa' (3.14) 

for a positive a'. Moreover, D*D < ^AHf and E*E < -^A 3 Hf (see, e.g., 
[GLL, Lemma A. 4]; note the different prefactors due to the missing ^ in 
their definition of ^4). Since A 2 > we get 

H P >\P- dT(k)\ 2 (l - 2aa) + H f (l- —A - -^-A 3 j - era' 

\ air 3a tt J 

> (\P\ 2 + \dT(k)\ 2 )(l-2aa)+H f (l-—A--^-A 3 -2\P\) - aa' , 

y air 6a tt J 

(3.15) 

where the last inequality holds for a < l/(2a). Since \P\ < \ by assumption, 
we can choose the constants a and a! appropriately and therefore arrive at 
(3.11). 

Using the results obtained above, we can now estimate (^q\Hp\^q), for 
an approximate ground state \l/o, from below. By Schwarz' inequality 

D 2 + D *2 > — aD* D - ^DD* (3.16) 

for any a > 0. Choosing a = l/s/a > \ we get, using DD* = D*D + ^A 2 , 

D 2 + D* 2 + 2D*D> 1 D*D-^A 2 . (3.17) 

Together with D*D < ^AHf and the a priori estimate {^ \Hf\^ } < 
const, a this allows us to conclude that 



(*oM 2 |tf > > " a2 " const - « 3/2 • ( 3 - 18 ) 



By similar arguments, using Schwarz' inequality as in (3.12), with a = 1/ yfa, 
and (v^ol \dT(k)\ 2 |* ) < const, a, 

\(V \y/adT(k) ■ A\W )\ < const. a 3/2 . (3.19) 
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We therefore obtain 

(*o|#p|^o> > \P\ 2 + -A 2 - const. a 3/2 + 

7T 

<*o| \<R{k)\ 2 - 2P • dT(k) + 2^P • A + yfea -B + H f + a 3/2 |^ > , 

(3.20) 

where we inserted the last term a 3 / 2 because it will be convenient later. The 
expression in the second line can be rewritten as 

£ [Wnl \dT(k)\ 2 + H f -2P- dT(k) + a 3 / 2 |Vn) 

n>l 

+ 2^K(ij n \2P-D* + (j-E*\ij n _ 1 )j , (3.21) 

where ip n denotes the n-photon part of the wave function ^o- We introduce 
the operator 

Q = \dT{k)\ 2 + H f - 2P • dT(fc) + o?l 2 , (3.22) 

which is positive for \P\ < ^. By using Schwarz' inequality, we can estimate 
(3.21) by 

(3.21) > -a(ip n -i (2P ■ D + a ■ E)^(2P ■ D* + a ■ E*) • (3.23) 



n>l 



We now investigate this term. The operator (2P '• D* + a ■ E*) acting on ipn-i 
creates a photon. Explicitly, in momentum representation, suppressing the 
polarization in the notation, 



1 n 

[£*Vn-i](fei,---fen) = —j= y^H(ki)tp n -i(ki,...,^i,...,k n ) , (3. 

v n 

v i=i 



24) 



where % means that k{ is omitted from consideration. An analogous expres- 
sion holds for D*. The right side of (3.23) can therefore be split into two 
parts, one coming from the "diagonal terms", where both on the left and 
on the right side the photon is created at the i'th position, and the mixed 
terms. More precisely, denoting F(k) = 2P • G(k)I C 2 + a • H(k) and using 
permutation symmetry, 



(3.23) = -a {In + Iln) , 



(3.25) 



n>l 



where the two terms I n and II n are given by 



F(k n )ip n -i(k 1 ,...,k n - 1 ) — F(k n )ip n - 1 {k 1 ,...,k n - 1 )j (3.26) 
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and 



II n = (n- 1)^(^)^-1(^1, • • • , kn-i^FfaWn-^h, fc„)) . (3.27) 



To estimate I n from above, we write Q = A + b, with A = \k n \ 2 + \k n \ — 
2P ■ k n , and use 

1 _ 1 b 
A + b A A(A + b) 1 ' ; 

and the estimates b > -2\k n \ J27=i A > \kn\0- ~ 2 \ p \) and „4 + 6 > 
\k n \(l — 2\P\) on the last term. Writing ip n -i = ip n -i{ki, ■ ■ ■ , k n -\) and 



W = W(h,...,k n ) = 



for short, we get 



2 1 k n | y] t -_i \ki\ 
\k n \ 2 + \k n \-2P-k n ' |fc n | 2 (l-2|P|) 2 



+ 



(3.29) 



In < (F(/c n )^ n _ 1 |vy|F(fc n )Vn-l 



= • G(k n )i> n ^\W P ■ G(fc n )Vn-i ; 

3 

+ J2( H j(kn)^n-i\w\H j (k n )^ n - 1 ) , (3.30) 

j=l 

where we used again the anti-commutation relations of the Pauli matrices 
and the fact that the mixed terms vanish, as explained in the upper bound. 
Inserting the expressions (2.8) and (2.9) for G(k) and H(k) yields (compare 
with (3.9)-(3.10)) 



_j_ / .3 4(|P| 2 |fc| 2 -|P-fc| 2 )+2|fc| 4 
(27r) 2 7| fc |< A |fc|3(|fc|2 + |A;|-2P-fc) 

4(|P| 2 |Af - \P ■ k\ 2 ) + 2\k\ 4 



In< ll^n-lf— U / <fk 

[2tt) z J\ k \<A 

+ (^-11^/1^-1)77^2/ d 3 k- 
{2tt) z J\ k \<A 



|fc| 4 (l-2|P|) 2 



. (3.31) 



To estimate the term II n we use the crude estimate 

li^Wn-ll < ( 2 I P • °( k n)\ + 3\H(kn)\) |Y>n-l(&l, . . . , fc„-l)| 

= \F{k n )\\ip n - 1 {k 1 ,...,k n - 1 )\ , (3.32) 

where the absolute values of tp n -i, G and H contain the appropriate norms 
in C 2 . On the n-photon space Q > (1 - 2|P|)(|fci| + \k n \ + a 3//2 ), so we get 
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the upper bound 



\lpn-l(ki, • • • , fc n _i)||F(fe n )||F(fci)||V'n-l(A;2, • • • , K 



(3.33) 



Introducing the photon density 

/ o n _i(fc) = (n-l) y a! 3 fc 2 ---d 3 A; n _i|V'„-i(A;,fe 2 ,...,fc n _i)| 2 (3.34) 
and using Schwarz' inequality twice, we estimate 

// n < 1 f d ^ K \/Pn-l{kl) P n-l{K)\F{k{)\\F{K)\ 



1-2|P|7 N + IM + a 3 / 2 

* T^p\^ h ^ [J d hd ""| fcl p„|(| fcl | + |M + « 3 / 2 ) 2 , 

(3.35) 

Since \F(k)\ < const, j/cl" 1 / 2 for small k the last term diverges logarithmi- 
cally as ol — ► 0. 

Putting together (3.31) and (3.35) and using the a priori knowledge that 
(*o|-f^/|*o) < const, a we get 

T(I+II)< J_f d3k H\PW-\p-kf) + 2\k\* 

|l { n + n) ~ (2vr) 2 y |fe| < A |*| 3 (|fc| 2 + |*| - 2P • k) 

+ 0(aln(l/a)) . (3.36) 
We insert this into (3.20) and arrive at the lower bound 

lril2 a 42 a /" 4(|P| 2 |fc| 2 - |P- £f) +2|fe| 4 

fp > \P\ + "A 2 - 7 T7T / C? 3 fc V ' ' 



(27T) 2 7| fc |< A \k\3(\k\2 + \k\-2P-k) 

- 0(a 3 / 2 ) . (3.37) 
Together with the upper bound (3.10) this proves (3.4), since 
1 1 



2\k\ 



\k\ 2 + \k\ —2P ■ k \k\ 2 + \k\ 



^ ( 2P . k+ I^L] , (3.38) 



(\k\ 2 + \k\) 2 V \k\ 2 + \k\-2P-k, 

and the first term on the right side vanishes after integration over k. □ 
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The coefficient of \P\ 2 in an expansion of Sp around P = determines 
the physical mass m as a function of mo, a and A. Since 



f 3 |fej /]Pj; 

4i<^ V + i*iUi 2 



:2m 



16vr 



In 1 + 



2 1 

(|£f + |fc|) 2 " jfcp 

A \ 3A(A+|m ) 



2m / 4 (A + 2m ) 2 



(3.39) 



Equation (3.4) leads to the following relation of m and mo, to leading order 
in a: 



m = mo 1 + a 



16 

3^ 



In 1 + 



A \ 3 A(A + | m 



2m J 4 (A + 2m ) 2 



(3.40) 



This relation is decisive when renormalizing the bare mass in the expression 
for the binding energy derived in the next section. Note that for any given 
positive values of m, a and A there is a unique positive solution mo to the 
equation above. This is in contrast to the analogous expression in dipole 
approximation, where the term in square brackets is replaced by a factor 
linear in A/mo, and consequently a positive solution only exists for a A small 
enough. More precisely, in dipole approximation the physical mass can be 
calculated explicitly, not only to leading order in a, and the result is 



m 



4a A 
mo + —A . 

37T 



(3.41) 



This was first shown by van Kampen in his thesis, as described in Kramers' 
biography [D, Sect. 16.III.D]. 



4 Binding energy 

In this section we shall calculate the binding energy of an electron in the 
field of nucleus of charge Z > 0. As we explained in the Introduction, we 
want to consider the case of small a with Za not necessarily small. To avoid 
confusion, we denote the fine structure constant appearing in the external 
potential by [3. I.e., we add to the Hamiltonian the potential 

V(x) = , (4.1) 

with P ~ 1/137. More general external potentials can be treated in the same 
manner, but for the sake of simplicity we shall not do so here. 
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Let eo = —^mo(/3Z) 2 be the ground state energy of the hydrogen atom 
without coupling to the photon field, and denote the corresponding ground 
state wave function by <j>Q. Let £(V) denote the ground state energy of (2.1). 
The existence of a ground state for H has recently been shown in [GLL]. 

The following Theorem gives £(V) to leading order in a. To obtain the 
binding energy, one has to subtract it from the self energy £(0), which will be 
done in Theorem 3. We introduce the notation that A(x) = D{x) + D*{x) 
and B(x) = E(x) + E*{x), where D and E denote the part of A and B, 
respectively, containing only annihilation operators. 

THEOREM 2 (Ground State Energy). With A denoting the operator 
A = 2^~\p\ 2 + V — e o + Hf, the ground state energy of (2.1) with external 
potential (4-1) is given by 



£(V > = e ° + ^ A2 -i^» ^JP^' + E^E- *,>+0(«*»). (4.2) 



Proof. We set again uiq = ^ for simplicity, and start with the upper bound: 
We use as a trial state 



* = <<h T, 



|p| 2 + y-e + |A;| 



2p-D* + a-E* 



6 T,0,... 



(4.3) 



with E*, D* and 4>o as explained above. Note that in contrast to (3.5) there 
is no need for an infrared cut-off, since the norm of ^ is finite. This can 
be seen as follows. An infrared problem can only arise from the overlap of 
p • D*(j>Q with 0o, since \p\ 2 + V — eo has a spectral gap above zero. Since 
p4>o is orthogonal to <fio, we can estimate the projection of p ■ D*(j>Q onto (po 
(for fixed k) by 



\{th\p-D*<h)\ = \G{k)-{<l>o\e- ik - x P <h 



= \G(k) ■ {<j>o\(e- ik - x - 1)W> )I < \G(k)\ \k\ ||W>olHl^o| 



(4.4) 



where we used that \e lk ' x — 1| < \k\\x\. The thus obtained additional factor 
| A; | makes the integral over k finite. 

Similarly to the calculation in the previous section, we obtain 



q( (2p ■ D* + a ■ E*)<j>Q T j (2p ■ D* + a ■ E*)<f> Q | ) + 0(a 2 ) . (4.5) 
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Using the anti-commutation relations for a, we see that the mixed terms 
vanish, and we arrive at (4.2) as an upper bound. More precisely, 

1 



and 



9 (Et<k 



A 
1 

A 



a ■ E*cf> T 







E*c 



, 



(4.6) 



(4.7) 



as can be seen by using the fact the G is real and H is purely imaginary, 
and that the operator A commutes with the reflection x — > — x. 

Lower bound: Proceeding as in the proof of the lower bound in Theorem 1 
(c.f. also [H]), it is easy to see that, in the ground state ^o> o\H f\^> o) < 
const, a. Moreover, 



v 3/2 



{^q\H\^q) > e + -A 2 - const. a a 

7T 

+ (* | \P? + V - e + 2^p ■ A + yfaa ■ B + if)|tf ) (4.J 



(compare with (3.20)). The expression in the second line can be rewritten 



as 



(V-ol \p\ 2 + V- eolV'o) + J2^n\ \p\ 2 + V-e + H f \^ n ) 



n>l 



+ 2^^{^ n \2p-D* + a-E*\^ 1 ) . (4.9) 



With A as stated in the theorem and F* = 2p ■ D* + o~ • E*, we get as a lower 
bound for (4.9), using Schwarz' inequality, 



(i/j \A\tpo) -aJ2(^ 



n>0 



F .±.r 



(4.10) 



To bound the terms with n > 1 from below, we need the following Lemma. 
LEMMA 1. Let N = dT(l) denote the number operator. Then 



F ■ —F* < const. (1+A0 . 



(4.11) 



Proof. By Schwarz' inequality 
1 



F ■ ^-F* < 2a ■ E—a ■ E* + 8D ■ p—p ■ D* . 

t/~\. t/~X t/~X 



(4.12) 



We first consider the last term. Let pi, i = 1,2, 3, denote the components of 
p. We claim that 



1 |eo| 

Pi , |2 it/ ~lT Vi - 1 + 77~ 

\p\ z + V-eo + Hf Hf 



(4.13) 
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It suffices to prove this for Hf some positive number fi. Since 

e P 1<e\p\ 2 < \ P \ 2 + V-e /(l-e) 
for < e < 1, we get 

1 1 

Pi \ p \2 + V-e /(l-e) Pi ~ e ' 



With e = [i/(fi + |eo|) this gives 
P 



i i 2 , T/ : — Pi < H 



(4.14) 
(4.15) 

(4.16) 



Using that A > Hf and again Schwarz' inequality for the terms with pi and 
Pj for i 7^ j, we thus obtain 



F • -F* < 2ct • E——a • F* + 24L> ■ 1 + 



.4 



F 



l£ol 



D* 



(4.17) 



Consider now the operator Di (^l + ^j-^j D* restricted to the (n — l)-photon 
sector. By use of Schwarz' inequality we can estimate 



D, ( i + M , D ; 



< n(V„-ie i/c ™' x Gi(fc n 



1 + 



IgoJ 



^e^G^kn)) , (4.18) 



where tp n -i is short for ip n -i(x, ki, . . . , k n -±). Since Hf > \k n \ in the 
n-photon sector, the last expression is bounded above by 



n 



! >n\\ 2 j2 I <Pk\G$(k)\ 2 (l + 

A=l J ^ 



|£o| 
|*| 



(4.19) 



Proceeding analogously for the other terms in (4.17) we prove the desired 
result. □ 



From [GLL, Thm. 6.1] we infer that in the ground state we have 
("fol-^j^o) < const, a, where the constant depends only on A. Using this 
we conclude from the Lemma above that 



(4.10) > {ifo\A\ifo) -aU 



const, a 



(4.20) 



and also that HV'oll 2 ^ 1 ~~ const, a. It therefore follows from first order 
perturbation theory [K] that, for small enough a, 



(4.20) > -aU j 



p. If* 
a 



t ) — const, a . 



(4.21) 
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The spin f is just one possible choice, any other spin will do. As already 
noted in the upper bound (c.f. (4.6) and (4.7)), 



Ait 



T) 



00 



4p ■ D- A p -D* +E- 
A A 



which finishes the proof of the lower bound. 



(4.22) 
□ 



With the expression for the ground state energy £ (V) in hand, we can 
subtract it from the self energy to obtain the binding energy. This generalizes 
a result in [H], where a lower bound on the binding energy was obtained. 
Moreover, the self energy was calculated, to leading order in a, in [H], and we 
will use the result from there. Note that both the ground state energy and 
the self energy are, to order a, linearly divergent in A (for fixed mo). The 
binding energy only shown a logarithmic divergence, which will be removed 
by mass renormalization in the next section. 

THEOREM 3 (Binding Energy). With B denoting the operator B = 
2^\p — k\ 2 + V — eo + | k\, the binding energy of (2.1) with external potential 
(4-1) is given by 



£(0) - £{V) = - e , + —Y J (P 



H<P0 



Gi 



+ — J V ( 00 ® Hi 

mi ^ 



i,j=l 

k-p 1 



k-p 



\k\ 2 + 2m \k\ B \k\ 2 + 2m \k\ 



b ®HiJ + O{a z/2 ) . 

(4.23) 



Proof. We set again mo = \ for simplicity. From [H, Thm. 1] we know that 



£(0) 



a 



A - a(0 



7T 

2a 



E(0) 



\k\ 2 + \kl 



E*(0) 0)+O(a 



= — (A-ln(l + A)) +0{a 2 ) , 



7T 



(4.24) 



where |0) denotes the vacuum in T^. We claim that the vacuum expectation 
value in the first line of (4.24) can be written as 



a( <pQ 



E- 



\p + k\> + V -e Q + \k\ 2 + \k\ 



E* 



(4.25) 



where E = E(x). To see this, consider the unitary transformation U on the 
one-photon space L 2 (M 3 , d 3 x; C 2 ) L 2 (R 3 , d 3 k; C 2 ) given by 



m){x,k) 



ik-x 



ip(x,k) , 



(4.26) 
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suppressing the dependence on spin and polarization in the notation. We 
have 

(UEt<f>o)(x, k) = Hi{k)Mx) =<t>Q®Hi. (4.27) 

Moreover, 

u(\p + k\ 2 + V(x) - e + \k\ 2 + \k\)u* = \p\ 2 + V{x) - e + \k\ 2 + \k\ 

= (|p| 2 + V - e ) ® 1 + 1 ® (|*f + |fc|) , (4.28) 

which proves the claim. 

Introducing the notation B' = \p\ 2 + V - eo + |A;| 2 + |fc| = B + 2p ■ k and 
using this unitary transformation, we have 



E ■ 



E 



|p + fc| 2 + ^-eo + |fc| 2 + |A;| 
1 1 



E* 



H 



B B' 

3 

E (^o ® ^ • *4 + ±2p ■ kUp ■ k± 



i=i 



(4.29) 



The first term in the last expression is zero by symmetry with respect to 
the reflection x — > —x. Moreover, B' Vo <8> Hi = (p ® #i/(N 2 + l&D- 
Hence, using Theorem 2 and (4.24), and the unitary transformation (4.26), 
we obtain, up to an error not bigger that order a 3 / 2 , 



3 1 
£ (0) - £(V) = -e + 4a ^ (p;0 o (g> G; - pjfo (g> G 



+ « E ® 



i=i 



2/c-p 1 2fc-p 



proving the Theorem. 



fli) , (4.30) 
□ 



5 Renormalization 



In the previous section, we calculated the binding energy to leading order in 
a. However, the bare mass enters this calculation, which is not a physical 
quantity. We have to replace it by the physical mass, which was, to leading 
order in a, calculated in Section 3. Note that both expressions show a 
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logarithmic dependence on A. It turns out the these divergences cancel and 
a finite result will be obtained. 

In the two terms in (4.23) of order a we can simply replace mo by m 
since we are only interested in a leading order calculation. This also affects 
eo and <fio; i.e., from now on eo = — \m{(3Z) 2 and 4>q is the ground state 
wave function of + V- 

The mo in the first term eo in (4.23) has to be replaced by the expression 
(3.40), however. This leads to 



m 
~2~ 



/ 16 


In ^ 









A\ 3A(A + |m) 



2m J 4 (A + 2m) 2 



+ — 9 > 

m z ^— ' 



i=l 



)ft 



k ■ p 



k • p 



\k\ 2 + 2m\k\ B \k\ 2 +2m\k\ 



</>0<8ft), (5.1 



with B = ^\p — k\ 2 + V — eo + \k\. We suppress higher order terms in a 
from now on. 

We shall show that (5.1) has a finite limit as A — ► oo, i.e., the cut-off can 
be removed. Let B' denote the operator B' = ^\p\ 2 + V — eo + ^\k\ 2 + \k\ = 
B + p ■ k/m. An easy calculation, using (2.8) and (2.9) and the fact that 



£G*(fc)G$(fc) 



(2vr) 2 |A:| 3 



(|fc| 2( % kikj} , 



(5.2) 



gives 

3 



^2 (pi<f>o ®Gi — pjcpo <8 G 



+ E 

i=i 



ft 



k • p 1 



|fc| 2 + 2m|A:| B' \k\ 2 + 2m\k\ 

'Mp\ 2 + V-eo A 



/ 



2m 



' 2m 



00 (8 ft 
Pj4>o 



(5.3) 



Here / is the function 

4 /" A 



A; 5 



1 

TT + 



1 



(5.4) 
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Note that Ylj {PjfiolPj&o) = 2m|eo|, and that 



/(0,A) = — ln(l + A)- 



3 A(A + |) 

4 (A + 1) 2 



Therefore (5.1) can be written as 



m 



£(V)-£(V) = ^Zf 



+ — y2(pj<f>o \f 



^\p\ 2 + V-e A 



j'=i 



2m 



^ "/ 0,— 



2m 



A \ 



2m J 



Pj<Po 



ij=i 

3 7 



i=l 



1 1 



\k\ 2 + 2m\k\ \B B'J \k\ 2 + 2m\k\ 



k • p 



(5.5) 



Hi 



(5.6) 



By definition B = B' — b with 6 = p • /c/m. However, it is easy to see 
that the expressions (5.1) and (5.6) do not change if we replace B by B' + b. 
Therefore we can replace l/B by \{\/{B' - b) + 1/(B' + b)), and instead of 
1/B — 1/B' in (5.6) we can write 

1 ( — X \ 1 - 1 1 1 f5 7l 



2\B'-b B' + bJ B' B> B'-bB'^b B 
For given operators h and p let Ck(h,p) denote the operator 

Ck ^ p) = hTWTW\ p - kx 

p-k- — — . (5.8) 



h + \k\ 2 + \k\ -p-k{h + \k\ 2 + \k\)~ l p ■ k h + \k\ 2 + \k\ 
We define the operator valued matrix T A by its components 



rMh, P ) 



2vr 2 L 



'|fc|<A 

being an operator on L 2 (R 3 , d 3 x; C 2 ) <8> C 3 . Using again (2.8), (2.9) and the 
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orthogonality relations of £\{k) and k, we can rewrite (5.6) as 



£(o) - £(v) = -(pzy 



*3 



2m 2m 
A/(2m) ( akbl 2 + - e p 



2m 



m 



Vj<t>o / 
Mo) • (5.10) 



Both terms on the right side have a nice limit as A — > oo. We define 



S(e)= lim /(0,A)-/(e,A) 



A-^oo 



4 r 

3vr io 



:fc 5 



{k 2 + k){e + k 2 + k) \k* (k 2 + k) 2 



1 

TT + 



and T(h,p) by 



T ij (h,p)=T*>(h,p) . 



(5.11) 



(5.12) 



Note that this operator, for h = ^\p\ 2 + V — eo, is well defined. To see this, 



we estimate 



2ml 

1 



< 



B'-bB'-H ~ \k\ 



(5.13) 



and 



JL(p . k f < \\p\ 2 \k\ 2 < \k\ 2 — {±\p\ 2 + V- eo/(l - e)) (5.14) 



for all < e < 1. Since B' > |£:| 2 /(2m) + |/c| we get 



T ( + L) < A /c „ (AW . + V- _ eo/(1 _ t)) " (5 . 15 ) 

\ 2m m / em v ^ m ' 9ir 



as an operator on L (R , dx; C )<8>C . The last factor comes from the integral 



d^k- 



\k\< 



16 

3ir 2 J K r '\\k\ 2 + \k\) 2 \\k\* ' (|fc| 2 + |fc|)V ' 



1 1 

+ 



(5.16) 



We see that, for h = +V — eo, both the operators T (h/2m,p/m) 

and /(0, A) — f(h/2m, A) are monotone increasing in A, and bounded by 
const. (1 + h) independent of A. Passing to the limit A — > oo, we have thus 
proved: 
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THEOREM 4 (Renormalized Binding Energy). To leading order in 
a, the renormalized binding energy, after removing the cut-off, is given by 



£(0)-£(V) = -WZf--J2{Mo 

£ III * 

3 

m \ 




PM , (5-17) 



with S andT the positive operators defined in (5.11), (5.12) and (5.9). 



Note that, by scaling, the right side of (5.17) is m times a function of 
(3Z. This is of course clear from the viewpoint of physical dimensions, since, 
after removing A, m is the only energy scale in the problem. 



The function 5 is monotone increasing and concave, with 

4 



and 



5(e) 



5(e) 



3vr 



eln(l/e) for e < 1 



h ln(e) 



for e>l 



(5.18) 



(5.19) 



Because of the logarithmic factor appearing in (5.18), the right side of (5.15) 
is, for eo <C m, much smaller than the second term in (5.17). I.e., for 
Z <C 1/(3 « 137, the term with T is negligible compared to the term with 5. 

Note that if we neglect T and approximate 5(e) by ^eln(l/e) we ob- 
tain exactly Bethe's result for the ground state energy shift [Be]. For this 
result it is not necessary to include the B field in the Hamiltonian, since 
the logarithmic factor in 5(e) for small e is entirely due to the term l/|fc| 4 
in (5.11), which comes from the p ■ A term in the Hamiltonian. The part 
with l/(|/c| 2 + \k\) 2 , which stems from the a ■ B term, is, when divided by e, 
bounded as e — > 0. 

For (3Z not too small, there is a significant difference between our ex- 
pression (5.17) and Bethe's formula. First of all, there is a different behavior 
of 5(e) for large e, and secondly the term T contributes, with a different 
sign than 5. It would be interesting to evaluate the terms in (5.17) numeri- 
cally, given the physical values of a = (3 and Z. For small Z, the result will 
essentially agree with Bethe's expression, whereas for larger values of Z our 
formula should be closer to the experimentally observed value. 



Remark. The appearance of the second term T in Theorem 4 is essentially 
due to the fact that the physical mass was obtained at total momentum 
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P = 0, whereas the ground state of an electron in the field of a nucleus 
shows a definite momentum distribution. In fact, a definition of the physical 
mass of the state 4>q by the relation 

(4> Q \£ P - £ o |0o> = 2^<to| \p\ 2 \4>o) , (5.20) 

where p = — iV x is the electron momentum operator, would account for a 
subtraction of a term similar to T. From a physical point of view, however, 
this definition is not very satisfactory since one wants to compare the binding 
energy with the rest mass of the electron and not with a mass defined in 
dependence on the state of the system. 



To see the order of magnitude of the shift of the binding energy due to 
the presence of the quantized radiation field to leading order in a, we now 
present a rough lower bound. The last term in (5.17) is positive and can be 
neglected. For the first term, we use concavity of S and Jensen's inequality 
to get 



m 



£(0)-£(V)>^((3Zf 



Q 

m 



3=1 



1 ELl (Pj<h 2^M 2 + V ~ e Pjfo 



2m 



E,=iWoMo) 



Using that Ei=i(Pj^o|Pj0o) = 2m|e | and 



(5.21) 



3 



1 



E 

j'=i 



Ph 



P 3 ^\p\ 2 + V-eo_ 
= 2vr(/3Z)|</» (0)| 2 = 2m 3 (pZ) 4 



we get 

In nature a = (5, and therefore 



8(0) - £(V) > ^([3Zf - ma(pzfS(((3Zf) . 



£(0) - £{V) > ^((3Zf [l - 2f3S{{f3Zf) 



(5.22) 
(5.23) 

(5-24) 



Inserting (5 = 1/137 and Z = 1 for the hydrogen atom gives a lower bound 
on the shift of the binding energy of —12.09 * 10 3 MHz, which is off the 
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experimental value by a factor of about 1.5 (see, e.g., the appendix of [Ki] 
for a detailed discussion). This is of course only a rough lower bound, the 
true value given by (5.17) is probably much closer. It should be noted, 
however, that the shift considered here is only the one due to the quantized 
radiation field, which is smaller than the shift due to relativistic effects. 



Due to the presence of the quantized photon field, the ground state is the 
only eigenstate of the Hamiltonian [BFS]. All the excited states that exist 
without radiation field turn into resonances. These are metastable states 
that decay after a characteristic lifetime. 

A reasonable guess for the energies of these resonances to leading order 
in the coupling a is to compute the expectation value of the Hamiltonian 
in a state similar to the one we used for the ground state and that has 
proved to give the right answer to leading order. I.e., we consider a state 
containing only one photon, but replace 0o an d eo by some excited state of 
the unperturbed hydrogen atom with corresponding energy. We denote by 
(fi n j an eigenfunction of 2^M 2 +^ with principal quantum number n, angular 
momentum Z, and corresponding energy c n (which is, of course, independent 
of I). Doing this calculation and also the mass renormalization to leading 
order in a, it is therefore natural to expect an energy of the metastable 
excited states, to leading order in a, as 



Note that the operator t^M 2 + ^ ~~ e n + 2ml^| 2 + \k\ i s n °t positive, but 
nevertheless invertible for almost every k. The expression (6.1) makes perfect 
sense if the integrals over k in (5.11) and (5.9) are interpreted as the Cauchy 
principal value. 

The same discussion as for the ground state energy shift in the pre- 
vious section also applies here. Neglecting T and approximating S(e) by 
^eln(l/e) gives the result of Bethe [Be]. This approximation is valid for 
Z <C 1/a. For larger Z, our formula presumably provides a better descrip- 
tion of the energy level shift due to the presence of the quantized radiation 
field than Bethe's formula. 



6 Lamb Shift 
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The correctness of the formula for the shift of the energy levels can best 
be tested on the classical Lamb shift, namely the splitting of the 2s±/ 2 and 
2p 1 / 2 states, since this is an effect entirely due to the quantized radiation 
field, and there is no splitting due to relativistic effects. 
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